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We formulate a simple model for a gas of extended hadrons at zero chemical potential by taking 
inspiration from the compressible bag model. We show that a crossover transition qualitatively 
similar to lattice QCD can be reproduced by such a system by including some appropriate additional 
dynamics. Under certain conditions, at high temperature, the system consist of a finite number of 
infinitely extended bags, which occupy the entire space. In this situation the system behaves as an 
ideal gas of quarks and gluons. 
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| I. INTRODUCTION 

o t 

The phase transition of strongly interacting matter has been intensively studied for many years. As early as 1960s, 
before the discovery of QCD, there was speculation about a possible new phase of strongly interacting matter, based on 
studies of the thermodynamics of a hadron gas. Particularly, in the Statistical Bootstrap Model the asymptotically 
exponential mass spectrum of hadrons implied the existence of a limiting temperature of about 170 MeV (the Hagedorn 
\ temperature) above which hadrons cannot exist. 
\Q • After the discovery of QCD and in particular asymptotic freedom, the discussion focused on the ideas of a Quark 
t-H \ Gluon Plasma (QGP), a system of weakly interacting quarks and gluons, and a possible transition between a pion 
gas and a quark gluon plasma. The physical picture for such a transition was that at the critical temperature the 
additional degrees of freedom carried by the quarks and gluons would be released leading to a rapid increase of the 
entropy, energy-density and pressure. This was supported by first Lattice QCD (LQCD) calculations 0. Lately, 
i— h ' however, the hadron resonance gas experienced a renaissance mainly due to its successful description of hadron yields 
in heavy ion and also elementary particle collisions [!, 0, H, H, 0, H, H, OH [HI Q3 • In addition it was realized that 
the sharp rise of the entropy density near the transition temperature, as observed in lattice QCD calculations, could 
be accounted for by a hadron resonance gas as well [13, LLJ]. The success of the hadron resonance gas below the 
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transition temperature T c , however, changes the physical interpretation for the QCD phase transition. Instead of 
releasing additional degrees of freedom at T c , the system has to reduce the number of active degrees of freedom at 
the transition temperature, because a hadron gas has many more than a Quark Gluon plasma. In case of a Hagedorn 
exponential mass spectrum, the entropy density actually diverges at the Hagedorn temperature 1 . But even if the 
number of hadrons is restricted to those needed for a successful description of hadronic final states in heavy ion and 
—i ■ elementary particle collisions, the entropy quickly exceeds that observed on the lattice. 

Meanwhile many papers have addressed this issue ranging from QCD-based approaches (see for example [III [III [III ) 
to various generalizations of the hadron gas [TtI EH O HO, HH, H2, [H, [24], [SfJ, l2ol [27| . In this article we will study 
the transition region, and we will provide a simple and intuitive modification to the hadron gas model in order to 
qualitatively reproduce the crossover transition observed in Lattice QCD [28| . Our calculations are based on the ideas 
proposed in [22j , where it was shown that under certain circumstances a gas of extended hadrons could produce phase 
. , transitions of the first or second order, and also a smooth crossover transition that might be qualitatively similar to 
^ ■ that of lattice QCD. 

We first observe that for the typical particle density at the transition the size of hadrons needs to be taken into 
account as it leads to a considerable suppression of the available phase space. As a result the number of effective 
degrees of freedom is reduced. This suppression alone, however, does not explain the ideal gas behavior of lattice QCD 
at high temperatures, above T c . Additional model assumptions have to be made. For example one could explicitly 
introduce a deconfined phase and then match the two different phases. For this exercise to work, however, rather 
detailed assumptions about the intrinsic nature of hadrons have to be made in order to avoid the occurrence of an 
actual phase transition [l^]. In this paper, we follow a different route. We adopt the philosophy that the same 
partition function should describe both "phases" . To this end we need to introduce appropriate dynamics in order 
to model the crossover observed on the lattice. This approach is similar in spirit to [22j. We find that the MIT bag 
model 29] of the hadrons is well suited for our purposes, because it embodies confined and deconfined phases from 
the very beginning. Thus, we will describe hadrons as extended bags of QGP and we will assume an infinite mass 
spectrum of the Hagedorn's type. The additional dynamics needed to describe the transition are simply the elastic 



1 Actually, for given choices of the mass spectrum parameters the entropy is finite at the Hagedorn temperature, but it diverges at any 
higher temperature. 
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interactions between hadrons. They give rise to a kinetic pressure that in turn "squeezes" the bag-like hadrons. We 
find that the behavior of our system depends sensitively on the choice of parameters for the mass spectrum. One can 
obtain either a real phase transition [2(1 or, as we shall show, a crossover. In addition, in the latter case the 
specific parametrization for the mass spectrum affects the microscopic structure of the gas of "compressible" hadrons 
at high temperatures (T > T c ). One finds either a high temperature phase that is populated by one or few infinitely 
large bags, consistent with the usual picture of a QGP. Or, for a different choice of parameters, one obtains a system of 
many, densely packed heavy hadrons 2 , which nonetheless exhibit the thermodynamic properties of a QGP of massless 
quarks and gluons. 

On first sight, our approach appears to be similar to the ideas of percolation models [3l|, HH, [H, [H, HH, [36|. The 
finite size corrections to the statistical ensemble remove all the configurations with overlapping hadrons, resulting 
in large hadronic states dominating the partition function. This is similar to percolation. However, there are quite 
some differences in the specific implementation. Fi rst, in our model, we do not consider an explicit coupling between 
the bags as it is done in the percolation model of [H, [H, HH, HI] . Instead, we take the effect of the kinetic pressure 
onto the bag sizes into account, resulting in a self-consistency relation for the effective bag pressure. This is more in 
the spirit of a mean field description, although we do not introduce an additional interaction but simply consider the 
kinetic pressure. Second, in contrast to purely geometric percolation, in our model the number of hadrons and their 
sizes are not independent quantities. In a given multihadron state, melting two or more hadrons together (to form a 
bigger one) results in a different kinetic pressure and, in turn, in the rearrangement of the sizes of all the hadrons. 

Throughout this paper, we will maintain a simple schematic approach to highlight the main features of the model 
leaving a more detailed quantitative analysis and further generalizations to future work. We will confine ourselves to 
the simplest case of nonrelativistic Boltzmann particles and we will neglect subtleties such as surface effects and van 
der Waals type residual interaction among the bags. 

This paper is organized as follows: in Section |TT] we introduce the main ideas of the model and derive the grand- 
canonical partition function for the gas of compressible hadrons. In Section [IIII we will use the corresponding isobaric 
partition function to perform a comprehensive numerical analysis. We will study the pressure, the energy density and 
the entropy density of the system. We will further analyze particles number, the filling fraction and the average mass 
of particles in the system. 



II. THE PARTITION FUNCTION 



In this section we will set up the general formalism for our model. Let us start with the partition function Z(V, T) 
of an ideal gas of Boltzmann particles of mass m and degeneracy g in the nonrelativistic limit. For the subsequent 
discussion it is advantageous to express the partition function Z(V,T) in a multiplicity expansion, i.e. as a sum of 
partition functions Z^(V,T) for fixed particle numbers N: 
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Here V and T are the volume and the temperature of the system, respectively. The function Z is the grand- canonical 
partition function with vanishing chemical potentials. In the context of this paper we shall refer to Z/v as the 
canonical partition function keeping in mind that this notation deviates from the conventions for relativistic hadron 
gases, where the canonical ensemble has fixed Abelian charges (such as electric charge, strangeness, baryon number), 
but no constraints on the number of particles. Eq. |T]) can be easily generalized to a multispecies gas of particles. If 
we label with j = 1, . . . , K the various particle species we have: 
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2 In what follows, we will sometimes use the word "hadron" with its widest mea ning without distinguishing among hadronic state such 
as resonances, bags or very short living states such as clusters (see, for example |30|0 . 
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In case of K — > oo, it is convenient to replace the discrete index j with a continuous spectrum density p(m) so that 
the number of species in the mass interval [to, m + dm] is given by p(m)dm. Formally, we make the substitution: 
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By expanding the exponential in Eq. ([3]) the partition function can then be written as: 
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Because a hadron gas, or, more precisely, a Hagedorn gas is characterized by an exponential mass spectrum, we set 
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where the parameters cq and a will be determined from empirical data. In the case of gas-of-bags models, which 
also have an exponential mass spectrum, the parameters cq and a will have to be determined from the underlying 
(dynamical) model parameters, such as bag pressure, and so on. We note that Co has dimensions of [mass]" -1 and a 
typically ranges from a = to a ~ 7 depending on the model. In Eq. ^ , Tq simply parametrizes the mass spectrum. 
In the context of the MIT bag model, To can be interpreted as the effective "temperature" inside the bag, as will be 
discussed in section |H"A1 By substituting Eq. © in Eq. §5$ we arrive at the following partition function for a hadron 
gas 
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This partition function Z(V,T) (and also Zn(V,T)) is divergent for T > To as already pointed out by Hagedorn. 
Although an upper limit in the mass spectrum p(m) regulates the divergences, it will not prevent the system from 
having a much higher entropy density than that observed on the lattice. 

An exponential mass spectrum without any cutoff may certainly be an oversimplification and a more realistic 
calculation may take into account discrete states as well as a mass spectrum that grows less than exponential above 
a certain mass. However, empirically the known hadronic states do indeed grow exponentially up to a mass of 
to ~ 2 GeV. Above that, very few states are known and it is not clear if this is an indication of a saturating density 
of states or simply the lack of experimental data on higher mass resonances. Therefore, working with an exponential 
mass spectrum without any cutoff appears to be an approximation as good as any other. Furthermore, because we 
are interested only in bulk thermodynamic quantities such as energy density and pressure, the use of a continuous 
mass spectrum should be a reasonable approximation as all these quantities represent integrals/sums over the mass 
spectrum. Therefore, in this paper we will assume that the mass spectrum is of the Hagedorn type and will discuss a 
dynamical scenario in the framework of the MIT bag model, which will regulate the partition function. 



A. The regularized partition function 



To develop the partition function of our model, we need to recall some of the basic features of the MIT bag 
model [29] ]. In its simplest formulation, hadrons can be considered as bags of partonic fields confined in a spatial 
region with a constant potential energy per unit volume B, where B is commonly referred as the bag constant or bag 
pressure. The total energy, i.e. the mass to of a bag with volume Vb is then given by [29] : 

m = U + BV b , (8) 

where U is the internal energy of the field inside the bag. When its linear extension is larger than the wavelengths of 
the partons (the quanta of the inner field) , we can approximate the bag by a gas of free massless particles confined to 
its volume [291 ] . For a sufficiently large \%, the internal energy is then given by the relation: 



U = 3p r V b , 
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where p r is the pressure of the gas. For a single hadron the stability condition requires: 
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which then gives 



m = ABV h 



(11) 



The effective temperature of the bag is related to the bag pressure by the relation To = kB 1 / 4 , where A; is a dimension- 
less constant whose value depends on the number of internal degrees of freedom of the gas inside the bag. For large 
Vb, the entropy S of a bag, is the entropy of a massless gas with internal energy U and pressure B [29j . therefore 3 , 
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From Eq. (|12p one can derive the level density 
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Note that the generic spectrum introduced in Eq. ([6]), has an additional contribution: mT a . This factor can be 
interpreted as a logarithmic correction to the entropy of the bag. In what follows, we will retain the spectrum in 
Eq. ^ and we will analyze different values of a. Of course, Eq. (fT3")) will correspond to the case a = 0. 

Once the temperature T of the gas of bags approaches To, the average masses and hence the volume (see Eq. (JTTJ)) 
of the bags grow very fast. Therefore, the bag-like hadrons tend to occupy more and more of the available space 
and, eventually, they will overlap. To avoid multiple counting of the phase space, configurations with overlapping 
bags need to be excluded from the partition function. For a finite system of volume V this can be achieved with an 
excluded volume correction, where the total volume V is replaced by the available volume (V — YliLi Vi) N > where Vt 
is the volume of the i—th particle. In addition the volume of all bags X)i=i ^ should not exceed the total volume V. 
Following Ref. [2(j, this leads to the modified iV-particle phase space integral 
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This modified phase space integral results in a well-defined and finite partition function at every temperature 
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It can be shown that such a system of extended hadrons leads to a constant value for the energy-density e, in contra- 
diction with Lattice QCD, where the energy density is found to increase with the fourth power of the temperature, 
^Lattice ~ T 4 . In fact, for T — > oo, the most favorite configurations are those where the hadrons occupy all the 
available space. In this case, the energy density of the system correspond to e = 4£>, i.e. the inner density of the 
hadrons (see Eq. (fTTjl ). The underlying reason for this behavior is that the system is not able to pick up additional 
kinetic energy once the entire volume is filled with bags, as the bags have no more room to move. 

Obviously some additional dynamics needs to be included to allow for the system to pick up more energy as the 
temperature is increased. To this end we adopt the idea of compressible bags |23|. More precisely, we will allow 
the volume of the hadrons to vary under the effect of the pressure generated by their own thermal motion in a self- 
consistent way. Consequently, as the temperature and hence the pressure increase, the bags will be compressed and 
acquire a higher internal mass/energy density. We will show that the system does not exhibit any limiting value of 
the energy density, and, under appropriate conditions, exhibits the desired increase of the energy density and entropy. 



3 In principle, on the left-hand side of Eq. 1 1 L'l ) one should subtract a constant So which corresponds to the entropy at U = 0. Here, this 
term has been omitted as it is immaterial for our purposes. 
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To illustrate the underlying mechanism, let us consider a gas of many hadrons. For small temperatures, T <C To, the 
system is dilute (V 3> J2i=i ^5) anc ^ behaves like a gas of noninteracting point-particles. With increasing temperature, 
the average mass, and hence the spatial extent of the hadrons, increases and as T approaches To the dilute-gas 
approximation seizes to be valid. The pressure exerted by the other particles becomes sizable and its effect on the 
hadrons properties, such as the size, can no longer be ignored. In other words, in addition to the bag pressure B, every 
particle in the system will feel an additional kinetic pressure pk that is generated by the thermal motion of the other 
hadrons in the gas. In this situation, the stability condition, Eq. (|10p . needs to be modified by taking onto account 
the contribution of the kinetic pressure pk ■ Microscopically, the pressure pk can be interpreted as the consequence of 
elastic collisions. The effect of inelastic collisions, which are certainly present in a hadron gas, in our approach are 
accounted by the infinite mass spectrum of hadrons without enforcing any constraint on the number of particles N. 
In this way all the possible configurations with few large hadrons or many small ones are included. This is analogous 
to the hadron-resonance gas model, where a large part of the inelastic interaction is taken into account by adding 
resonances as free particles in the gas. 

Neglecting any surface effect, the simplest generalization of the stability condition, Eq. (fT0|) . is 

p r = B+p k {V,T) . (16) 
A pictorial illustration of the pressure balance in the last equation is given in Fig. ([T]). Each hadron in the gas is 




FIG. 1: Pictorial representation of the pressure balance in a gas of compressible hadrons. The internal hadron pressure must 
be equal to the sum of the constant bag pressure B plus the pressure pk generated by the thermal motion of the hadrons 
themselves. 



characterized by the same internal pressure p r . But instead of being a constant (as in the case of a single hadron in 
the vacuum where p r = B), p r now depends on T and V , and must be evaluated in a self-consistent fashion from the 
partition function itself. 

It is clear that the number of hadrons, their sizes and the kinetic pressure are all connected by the above self- 
consistency relation. Accordingly, as we have already pointed out, if we split or combine two or more hadrons, the 
pressure, and thus the volume of all the hadrons in the system, changes. This new state corresponds to a distinct state 
of the ensemble that cannot be obtained by a simple geometrical clustering procedure as usually done in percolation 
models [H M, M, M, El . 

To account for this additional dynamics, we need to generalize the partition function, Eq. (|15|) . To this end, we 
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write the explicit dependence of volumes and the masses of the bags on the pressure p r 
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leading to the substitution 
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Finally, because the bag masses depend on p r , instead of integrating on dmi . 
the integral over the internal energies, i.e. we substitute 
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The factor 4/3 in the previous formula ensures that we recover the partition function, Eq. I|15j) in the limit of pk — > 0, 
i.e. in the dilute gas limit. With the replacements in Eq. (fT7|) to (|20|) . the modified grand-canonical partition function 
can now be written on the basis of Eq. (JTSJ) and reads: 
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The new partition function Z(V,T) is identical (by construction) to Eq. (115j) when p r — B. Conversely, for a given 
set {Ui}, the effect of a finite kinetic pressure pk > is to squeeze each bag to a smaller size (see Eq. (|T7)) ). and as a 
result, the effective bag temperature Tf, = kpl^ 4 increases. The Eq. (|2ip can be made more familiar by substituting 
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Obviously, in the dilute gas limit, pk — > % — > m,. 

We further introduce a lower bound m c for the integrals over {d?7i}. This is needed because, for a > 0, the spectrum 
in Eq. ([6]) has a pole in m = 0, resulting in a divergent partition function for a > 5/2. Because there are no hadrons 
lighter than the pion, we will set m c = m, = 0.139 GeV. 
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III. THE ISOBARIC PARTITION FUNCTION 



Because the pressure pk is thermally generated, it must be calculated from the partition function itself, resulting 
in a self-consistency relation. This is best achieved by introducing the isobaric partition function, which is defined as 
the Laplace transform of Z(V, T) over the variable V: 



Z(T,s)= dVZ(V,T)exp[-sV] 
Jo 



(24) 



The quantity sT in Eq. (|24[) plays the role of a constant external pressure. Accordingly, the equilibrium condition 
requires 



Pk = sT . 

The integral in Eq. (|2~4")1 can be solved analytically (see Appendix A) and gives: 
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In the limit V — > oo the asymptotic behavior of Z(V,T) is defined by the singularity of Z(T,s) with the largest real 
part [20(. We have two distinct cases: a < 5/2 and a > 5/2. For a < 5/2 the transform Z(T, s) in Eq. (f2"6"| has two 
kinds of singularities: the first, sq(T), is given by the pole of l/(s — /(T, s)), 



i.e. 



s (T) = /(T, s (T)), 



(28) 



whose solution is the pressure of the system according to Eq. (|25|) . The second singularity, s/(T), corresponds to a 
divergence of the function / itself. This happens if the exponent of the integrand of f(T, s) in Eq. (f27f vanishes, i.e. 
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The situation is schematically represented by the leftmost curve in Fig. ([2|). The solid line represents the function 
f(T, s) for a given temperature and the 45° 's dashed line corresponds to s. The intersection between the dashed and 
the solid line corresponds to the solution s = sq(T) of Eq. (|2"B")) and is denoted by a black dot. The function f(T, s), 
is a positive function of s that goes to infinity for s — Sf(T) and tends to zero as s — > oo. Consequently, there is 
always a solution for Eq. ([55]) and the pole s = s (T) = p k (oo,T)/T corresponds to the rightmost singularity. Thus 
the pressure has always a solution. For a > 5/2 the situation is different, however. In this case the function / has 
an essential discontinuity at Sf(T): it is finite at s — Sf(T), and diverges for s < s/(T). Since Sf(T) increases with 
temperature (see Eq. (|2T))0 and f(T, s) — > as s — > oo, for sufficiently large T, f(T, s/(T)) < s/(T), and, consequently, 
Eq. (f28|) does not have a solution. This situation is illustrated by the rightmost curve in Fig. ([2]) where f(T,sf(T)) 
(denoted by a X) lies below the diagonal 4 . Both these cases, have been discussed in [2(| [22j where the absence of the 
solution sq(T) was identified with the onset of a phase transition. Here, we analyze in detail the case of a crossover 
transition, i.e. a < 5/2. 

Before we proceed, let us fix the model parameters. In what follows, we set k = 0.68 and we keep the product 



T = fcB 1/4 = 0.17 GeV 



(30) 



fixed. This yields a bag pressure B = 3.9 • 10" 3 GeV 4 (i.e. B 1 / 4 = 250 MeV) which is a plausible value for the bag 
model parameter. The constant k has been chosen in order to obtain roughly the same value for e/T 4 as LQCD for 



4 Notice that this can only occur for temperatures T > kB 1 ' i = To since the singularity Sf(T) < for T < To (see Eq. 11291 ). 
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FIG. 2: Schematic representation of two possible curves f(T,s). When a < 5/2 the solution of Eq. (128f) always corresponds 
to the rightmost singularity (denoted with the black dot), whereas for a > 5/2 the function f(T,Sf(T)) is finite, and when 
T > To it can happen that s/(T) is the only singularity (denoted with the X). The figure is adapted from [2p|. |22||. 



a 





1/2 


1 


3/2 


2 


5/2 


c (GeV " 1 ) 


0.157 


0.199 


0.252 


0.318 


0.400 


0.502 



TABLE I: Standard values of the parameter Co for different choices of a. 



large T (see Fig. ([3])). Its value can also be estimated by counting the degrees of freedom of a thermal system of 
independent quarks and gluons . Our choice lies between the values for the lightest quark doublet u, d (k = 0.70) and 
for u, d and s quark (k — 0.66). The values for the remaining parameter Cq have been fixed by fitting the shape of 
the actual hadron mass spectrum over the mass range of 1 — 2 GeV. They are given in Table HI Of course a different 
mass range or a different choice for Tq would affect these fits. However, for our schematic considerations here, a fine 
tuning of the model parameters is rather meaningless. In the same spirit we also ignore the ~ 10% deviation of the 
LQCD result for e/T 4 from the free gas (Stefan-Boltzmann) limit [l4|. 

In Fig. (J3h') we plot the pressure Pfc(oo, T)/T A = so(T)/T 3 evaluated with the isobaric partition function, Eq. (|28[) . 
As one can see, the results depend on the choice of a. For a = and 1/2, the curves grow with the temperature with 
larger slopes for smaller a's. Instead, for 1 < a < 5/2, they settle onto constant asymptotic values (as we have verified 
numerically up to T ~ 6 GeV). As a change from 1 to 5/2 the asymptotic value converges very fast to 1/fc 4 = 4.67 
(the Stefan-Boltzmann limit 5 ) from above. As shown in the plot, the curves practically coincide with the Stefan- 
Boltzmann limit already for a = 2. This behavior can be understood by inspecting the solution Sq{T) = pk(oo,T)/T 
for the pressure. Because sq(T) is always larger than s/(T), we have 

p fc (oo,T)>^-B (31) 



5 In the situation where the system is completely filled by the inner hadrons matter (i.e. the free massless gas), the relation between 
pressure and temperature is p = (T/fc) 4 . Accordingly, the energy density is given by e = 3p = 3(T/fc) 4 . 
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FIG. 3: (Color online) Left panel (a): The ratio pk(oo,T)/T 4 calculated with the isobaric partition function for various values 
of a. Above To, the curves at a = 2.0 and a — 2.5 are practically coincident. Right panel (b): The corresponding value of 
e/T 4 . 



and for T — > oo 

Pk(<x,T) > _1_ 
T 4 ~ fc 4 



(32) 



The pressure pk(oo, T) converges to T 4 /fc 4 only for sufficiently large values of a, when the solutions sq(T) and s/(T) 
get closer and closer as T — > oo 6 . 

In Fig. (|3Jd), we plot the ratio e/T 4 , where e has been evaluated numerically by using the relation: 

£ = T w~ Pk - (33) 

Again this quantity converges to a finite asymptotic limit only for 1 < a < 5/2 and, as before, coincides with the 
massless gas limit 3/fc 4 ~ 14 for a = 2 and a — 5/2. The overall behavior is roughly the same as LQCD except for 
a small "horn" right above T . A closer look reveals that this is due to the constant bag pressure B that produces 
a contribution ~ B to the energy density of the system. Being a constant term its contribution to e/T 4 becomes 
negligible at high temperatures. 

In Fig. ([4^), the ratio (e — 3pfc)/T 4 is plotted up to T = 0.6 GeV. This quantity corresponds to the trace of 
the energy-momentum tensor O w (T)/T 4 , which is actually the fundamental quantity calculated in LQCD [l4j]. In 
Fig. iJUb); we plot the ratio s/T 3 = (e +pk)/T 4 , where s is the entropy density. Again, for a — and 1/2, the curves 
do not converge to a constant value, as expected from the previous results for e and pk- 

In our scheme, the pressure p^ (and thus the energy density and the entropy) exceeds the corresponding Stefan- 
Boltzmann limit, except for a = 2 and 5/2 where it converges to it. To obtain a lower pressure (at least at finite 
temperatures) one needs to further reduce the effective degrees of freedom of the system. This could be possibly 
achieved by introducing a surface energy term. In addition to the fact that part of the energy of the system would 
be spent to create the bags surface, such a contribution would favor spherically shaped bags, resulting in a further 
suppression of the accessible phase space. Probably, a similar effect could be also obtained including some residual 
repulsive interaction of the van der Waals type. 



Actually, since pk = so (T)T, to obtain the pressure of an ideal gas the difference (so(T) — Sf(T)) must decrease faster than 1/T. For 
the ratio p^/T 4 it is sufficient that the difference (sq(T) — Sf(T)) grows slower than T 3 . 



10 




We stress that the general behavior of our model for the pressure, energy- and entropy-density cannot be obtained 
by simply introducing an upper mass-cutoff on the exponential spectrum in Eq. |(7J): the entropy-density would 
considerably exceed that obtained from LQCD even if we kept only masses up to 2 GeV. In addition, the system 
would reach the massless gas limit only at temperatures much higher than the cutoff itself. On the other hand, in our 
model it is absolutely essential to assume an infinite mass spectrum. Otherwise the flat behavior in Fig. ([3]) and on 
Fig. (Hb) would be spoiled and all these quantities would decrease with the temperature. 

It is also interesting to plot the average bag effective "temperature" = k(p k + B) 1 / 4 . As shown in Fig. ([5]), for 
a = 2 and a = 5/2 this converges to T very quickly above To- For a = 3/2 or smaller the effective bag temperature is 
always larger than the system temperature. This fact is a direct consequence of the inequality in Eq. (|31[) that gives 

T b = k( Pk + B) 1 ' 4 >T, (34) 

Notice also that in the region T < T , T& ~ T for any value of a. In other words, the pressure p k is negligible 
with respect to B and the system behaves as a standard hadron gas. It is worth mentioning that in our framework 
the compressible hadrons can exchange energy only through a mechanical work (compression). They are, therefore, 
thermally insulated from the rest of the system and the bag temperature, as well as the entropy, must be understood 
as quantities that measure the degeneracy of the hadronic states. 

Finally, for 1 < a < 5/2, our model seems to produce a smooth crossover transition toward a new regime whose 
features are very similar to those of a gas of massless particles, even though no deconfincd states are included in 
the partition function. To better understand this behavior it is useful to study the particles density (n) = (N)/V 
(Fig. (j6ji)) and the filling fraction (/./.) (Fig. (JB]b)) which is defined as: 

i> r (^hadrons) 

/■/■ = y , (35) 

where (Vh a drons) is the average volume occupied by the hadrons (for a rigorous definition and formula see Appendix 
B) . The particles density can be calculated from the isobaric partition function by introducing a fictitious fugacity A 
(to be set to 1 afterward) for each particle in the system, i.e. replacing Z(V, T) with 

oo 

Z(V 7 T,X)=J2 * N Z N (V, T) . (36) 
Accordingly, the isobaric partition function in Eq. (121)1) becomes 
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FIG. 5: (Color online) The average effective bag temperature. The 45° 's straight line corresponds to the system temperature 
T. 



and the corresponding solution for the pressure pu 
can be then obtained as: 



Pk{oo,T, A). In the infinite volume limit, the particles density 



1 d\nZ(V,T,X) 
(n) = urn 

v 1 v^ooV dX 



£ dp k (oo,T,X) 
T ' 



dX 



(38) 



where in the last equality we have used the known relation InZ = pV/T. As shown in Fig. (|6^), the particles 
density grows very rapidly as T approaches To from below. This is qualitatively what one expects for the hadron gas, 
where the average number of particles shows a monotonically growing behavior. Conversely, for T > To this quantity 
depends very strongly on the choice of the parameter a. For a — 0, 1/2 and 1 we observe a change in the slope, but 
the curves still grow monotonically. For a — 1.5 and 2 the particles density has a local minimum at T ~ 0.2 GeV 
and T ~ 0.27 GeV, respectively (the latter lies outside the plotted region), and then it starts growing again with 
smaller slopes for larger values of a. For a = 5/2 the local minimum has disappeared, and after a sharp maximum at 
T ~ 0.18 GeV the particles density goes to zero. This is the effect of the finite size of hadrons, which tends to saturate 
the available system volume. In fact, as shown in Fig. (J6Jd), when T < To the filling fraction, Eq. (|35[> . is relatively 
small, whereas for higher temperatures, the system is almost totally filled by extended particles, i.e. (Vhadrons) ~ V ■ 
In this scenario, the space and the phase space available is strongly suppressed, and the system tends to be populated 
by a smaller number of heavy particles. This effect is strongest for a — 5/2 (the filling fraction converges very fast 
to 1) and therefore (n) — > 0. For smaller values of a this saturation effect becomes slightly less pronounced. A closer 
inspection reveals that for a < 1 the filling fraction has a maximum at very high temperature and then decreases with 
a very small slope 7 , whereas for a = 1 it seems to settle to a constant value /./. ~ 0.98. The phase space, therefore, 
is never entirely suppressed. 



7 For a = 0, the rate of decrease of the filling fraction for large temperatures is maximum, but still only ~ 2% going from T = 1 GeV to 
T = 10 GeV 
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The behavior of the system changes continuously by varying a. A numerical analysis indicates that there exist a 
value cto between 2.12 and 2.13 such that, at high T, the particles density vanishes for any ao < a < 2.5. In this range, 
the system is populated by one or few infinitely extended hadrons that occupy the entire space, filling the system with 
their inner QGP matter, which is a possible scenario for the deconfincd phase. Conversely, for 1 < a < ao, we find 
many, rather heavy, "squeezed" hadrons, which nonetheless mimic an ideal gas of massless particles. The number of 
particles, however, might be affected by the introduction of a surface energy term in the spectrum in Eq. Such a 
term would result in an energy cost associated with the splitting of a large hadron into many smaller ones and might 
then widen the range of values for a that lead to (n) ~ at temperatures above To. 

Note that, even though (n) can have a minimum or even vanish, the entropy of the system always increases 
monotonically with the temperature (see Fig. ([Us)). This is due to the fact that the dominant contribution (already 
at T ~ To) comes from the bags entropies. This has been checked by using the classical expression for a gas with (N) 
particles (with the excluded volume correction) and degeneracy g: 



S ~ S class {{N),g) = S class ({N),g = 1) + (N) Ing 



The contribution (N) Ing corresponds to the sum of the bag entropies in Eq. (p~8|) . and it is dominant at high T. 
Another interesting quantity is the average hadron mass (m), shown in Fig. ([7]) that has been evaluated according 

to 

H = A - §T (39) 
(n) 2 

As one can see, smaller values of a correspond to higher (m) at low temperature, whereas the situation is reversed at 
high T . This follows from Eq. (|3T)|) as a direct consequence of the behavior of (n) (Fig. (jS^)). For a — 5/2 the average 
mass (m) grows with the temperature with the maximum slope and diverges when (n) — » 0, i.e. when the system is 
populated by a finite number of infinite "hadrons" . 

In the vicinity of the transition region (to) varies from 4.34 GeV (for a = 0) to 0.75 GeV (for a — 5/2) at 
T = 0.16 GeV, and from 5.70 GeV to 1.08 at T = 0.17 GeV, for a = and 5/2, respectively. These values will 
somewhat depend on the choice of the model parameters. They would depend even more on an eventual upper mass 
cutoff. In fact, the high value of (m) for a = results from our assumption of an infinite mass spectrum and already 
at T — 0.16 GeV it falls outside the region of the known hadrons. Our estimates are, however, lower than what one 
obtains for an exponential spectrum of point-like hadrons as in Eq. ([7])- In such a case, for a = the average mass is 
(m) = 6.80 GeV at T = 0.16 GeV and infinity at T = T = 0.17 GeV (the partition function itself is divergent). A 
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FIG. 7: (Color online) Average hadrons mass. 



mass cutoff at 2.0 GeV would reduce these numbers to (m) — 1.36 GeV and (to) = 1.43 GeV at T 
T = 0.17 GeV, respectively. 



0.16 GeV and 



A. Consistency check 



As a final remark, we want to discuss the consistency of the model. To check this point we must make sure that 
quantities such as the energy density (that we have evaluated from the pressure by using Eq. ([33]) ) correspond to a 
thermal average of the form 



1 Estates ^state6 E " t "- t <>/ T 
v Z^states c 



(40) 



A first hint in this direction, is given, a posteriori, by the results shown in this section, particularly, by quantities such 
as the filling fraction. The /./. has been evaluated by making use of the relation in Eq. ([ST)]) (sec Appendix B) that 
implicitly relies on a form like Eq. (|40[) for the energy density. Indeed, a wrong thermodynamical interpretation of e 
would likely have lead to dramatic consequences on the filling fraction, which contrarily assumes only physical values 
in the interval [0,1]. However, to make a more direct test, we will provide an approximate expression for the grand- 
canonical partition function and we will compare e in Fig. (JSJa) with its corresponding value obtained as in Eq. (I40|) . 
This can be done for the case a — 3/2. For this particular value of a the multiple integrals over {drji} in Eq. ((23|) can 
be reduced to an unidimensional integral (see Appendix C) that greatly facilitates the numerical treatment. Although, 
as demonstrated in this section, results do depend on the choice of a, for the range ofO < a < 5/2, one can hope that 
the following arguments will be still valid. 

Using a = 3/2 in Eq. (|2"3"|) we get for the partition function 



Z(V,T)= J2Z N (V,T) 



(41) 



N=0 
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with 



Z N (V,T) 



x exp 




(42) 



Here we have added the suffix N to p r to indicate the dependence of the pressure p r ,N on the particle number in the 
canonical ensemble. The pressure p r> N is now just a parameter of the model. To determine its value, we need to find 
the value of p r> N that maximizes the logarithm of the integrand in Eq. (|42D : 



$ = 



1 



kpJ N 



- - 

T \4 



D 



^Pr,N 



N 



N 



Vi 



4:Pr,N 



(43) 



where (for the moment) we omit the function. The above expression depends on ^# = X)t=i Vi- To simplify the 
following derivation, we will introduce an approximation. Instead of solving d^/dp r _N — for a generic set rji, . . . , rj^, 
we find a solution for the most important configurations defined by the value ^ that maximizes the integrand. In 
other words, we solve the system of equations: 











(44) 



Note that the second condition ensures that e = T 2 (d\nZ/dT)/V has the form of Eq. (|40|) as any implicit dependence 
on T of Pr 7 N do not contribute to 9 In Z/dT. In fact, if we denote with d* \nZ/dT the derivative performed only on 
the explicit T dependence of In Z we have 



dlnZ d*lnZ 



dT 



dT 



dZ N dp TtN 



Z ' dv r N dT 

N=0 r ' 



where the second term vanishes because of the second condition in Eq. (|44p . The Eq. 
D): 



PrM = B 



NT 



V 



4p r ,JV 



(45) 

results in (see Appendix 
(46) 



The interpretation of the above expression is straightforward. Writing p r ^ — B + pk,N one obtains the equivalent 
equation: 



NT 



Pk,N 



V 



(47) 



4(B+p fc ,iv) 



Here, the right-hand side of Eq. (|47|1 is simply the canonical pressure of an ideal gas of N particles with the total 
volume V replaced by the available volume (V — Vi). It is then natural to identify pk,N with the kinetic pressure 
in the canonical ensemble. Of course, once averaged over N and over {r/i}, for a sufficiently large V , this pressure 
must coincide with p^{oo, T) evaluated with the isobaric partition function. Eq. (|47[) has the form of a self-consistency 
relation, as pk,N appears also in the right-hand side in the excluded volume term. Eq. (|47p is a quadratic form and 
has two solutions: a negative and positive one, and the negative corresponds to the situation where volume of hadrons 
exceeds the total volume, ^Zi=i Vi > V- The positive solution on the other hand ensures X)«=i V < V and therefore 

is always fulfilled. In what follows, we will adopt the positive solution of 

//x and we will integrate numerically 



the condition for the function in Eq. 
the Eq. (|47p for any {771, . . . , (not only for the most probable set 77 
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Zn(V, T) in the variables {r/i}. The partition function Z(V, T) is then evaluated by summing over the particle number 
N up to a cutoff N cut 

Ncv.t 

Z(V,T)=J2 z n(V,T) (48) 

JV=0 

where N cut is sufficiently large to ensure the accuracy of our calculations. Finally we test the consistency of our 
picture by comparing the ratios Pk/T and e/T 4 with the results obtained with the isobaric partition function. In 
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FIG. 8: (Color online) Left panel (a): Comparison between the ratio Pk/T 4 evaluated with the isobaric partition function (solid 
line) and with the grand canonical partition function (dashed line). Right panel (b): The same as in the left panel for the 
corresponding ratio e/T 4 . 



Fig. ©, the ratios Pk/T 4 (left panel) and e/T 4 (right) have been evaluated with the isobaric partition function (solid 
line) and with the grand canonical partition function (dashed line) for a volume V = 6.4 ■ 10 4 GeV~ 3 , which is, as we 
checked, a good approximation of the infinite volume limit. As one can see, in both cases the two curves are very close, 
actually for the energy density they are practically coincident. The small (expected) difference between the isobaric 
and the grand-canonical result reflects the quality of our approximation. The same test has been also performed on 
the entropy density, the particles density and the filling fraction. For all these quantities we have observed an equally 
good agreement. 



IV. CONCLUSIONS AND DISCUSSION 



In this article, we have studied the crossover transition of the gas of bags [22|. We have found that the behavior of 
the system depends sensitively on the parameter a of the Hagedorn-like mass spectrum p(m) = c^ra~ a exp {to/To}. 
The system exhibits a crossover transition for < a < 5/2 and an actual phase transition for larger values. In 
the range < a < 5/2 we made a coarse scan of a, setting a = 0, 1/2, 1, 3/2, 2 and 5/2. For 1 < a < 5/2 the 
gas of bags undergoes a sharp (yet, continuous) transition qualitatively similar to lattice QCD. In this range, the 
asymptotic values of p/T 4 , e/T 4 , s/T 3 , coincide with the Stefan Boltzmann limit for a — 2 and a — 5/2, whereas 
they settle to slightly larger values for a = 3/2 and 1. For a — and 1/2, these quantities grow indefinitely with 
the temperature with small, decreasing, slopes going from ev = 0toa = l/2. We have also studied the (strong) 
dependence of the particles density (n) = (N)/V (where V — > oo) on a. We have found that there exist a limiting 
value ao between 2.12 and 2.13 such that for ao < a < 5/2 the particles density vanishes at high temperature. The 
system is then populated by one (or few) infinite bag(s) that occupies the entire volume. Conversely, for a < ao, 
(n) grows with the temperature. In the range 1 < a < ao, the ideal gas behavior is mimicked by a number of heavy 
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extended bags that saturates the phase space forming a dense system. A pictorial representation summarizing the 
various high-temperature phases of the model is given in Fig. 



No ideal gas behavior 




Ideal gas behavior 
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2 *„ 
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FIG. 9: Pictorial representation of the various high-temperature phases of the model as a function of a. The scale has been 
distorted for visual reasons. 



In this work we have explored a simple, intuitive model for a gas of hadrons in the vicinity of the transition at 
vanishing baryochemical potential. To this end, we have adopted the idea of the MIT bag model and we have described 
hadrons as extended QGP bags. We have shown that, in the vicinity of To (which is directly related to the transition 
temperature of our model), elastic interactions among hadrons play a fundamental role. In our schematic approach, 
they are quantified by the thermal pressure pk- The effect of pk is to squeeze the hadrons, and for a certain set of 
model parameters, 1 < a < 5/2, the ideal gas behavior at high temperature can be reproduced. At the same time, 
the effective inner "temperature", or rather degeneracy parameter, T,, of the bags increases with pk, resulting in a 
temperature-dependent mass spectrum. Above To, the physical picture of the QGP phase corresponds to a number 
~ 1 of infinite bags that occupy the entire space. This is indeed the situation for ao < a < 5/2. On the other hand, 
a large number of independent QGP bags (such as for a < a$) would contradict the lattice findings of vanishing 
flavor-flavor correlations [13, [H, [3^]. However, in order to define precisely the range of values of a that lead to a 
consistent QGP scenario, it is fundamental to study the effect of a surface energy term. This contribution disfavors 
configurations with a large number of particles (they are more "expensive" in terms of surface energy) and might 
reduce the number of bags in the high-temperature phase, widening the range of possible values for a. 

Future work will likely concentrate on the effect of surface energy terms in addition to the study of van der Waals 
type residual interactions. It would be also interesting to analyze in detail the behavior of the system at finite 
baryochemical potential. This could be done starting from the formalism developed in [24| . 
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Appendix A 

By using the equations (| 16[) and (|25|) . the isobaric partition function in Eq. (|24[) can be conveniently written as: 

^ T > s )=EAnIT/ dr?i Mfo}' T '*) (49) 

N=0 Yi=l Jm ' 

x [°° dvexp[- s v] (v-Y—^ — ^ e(v-Y— 

L 1 \ fr{ 4(5 + sT) J \ Z^^B + sT)) 
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where 



, rp x 3A72 

h N {{ni},T,s) ee ( — J c, 



A' 



z— 1 v 



3/2-a' 



X exp ■ 



1 /3 



_k{B + sT) 1 / i T \4 A{B + sT) 
The integral over dT^ can be carried out 



N 



,oo f N 

/ dV exp [sV] I V - V] 



AT 



AT 



^ 4(5 + sT) 



^4(5 + sT)^ 



AT! 



A 



,AT+ 



T ex Pi-E 



- ll A(B + sT) 



and Eq. becomes: 



oo r Af „ 

%»4e n/ 



d^ 



feiv({7?i},r,5) 



N 



exp < - ^2 n 



N=0 Li=l ' 

The multiple integral in the last equation, can be factorized as: 



^ n 4(B + sT) 



JJy h N ({r)i}, T, s) exp | - ^ ^ ^ ^^ |e f(T, s) 



N 



One then obtains Eq. ||25J with /(T, s) given by Eq. (j2"7) . 



(50) 



(51) 



(52) 



(53) 



Appendix B 

By using Eq. (jTTJ) , we write the average volume occupied by hadrons as 

Madrons) = (N) SD. = {N ) ^ = (N) ^ , 

where we used the stability condition p r = p k + B. In turn, the average mass (m) reads (see also Eq. (|39p ) 

s 3 

The filling fraction /./. can then be expressed in terms of e, p k and the particles density (n) as: 

(^hadrons) £ ~ 3/2<n)T 



V 



(54) 



(55) 



(56) 



Appendix C 

The partition function of N particles can be reduced to a unidimensional integral when a = 3/2 as the factor 
(3/Arji + BVi) 3 ^ 2 a reduces to 1. Let us begin by making the substitution: r)i — > yf + m c that gives for Zn in 
Eq. 03): 



x exp ■ 



N 



1 

ivT 



1 1/3 B 



TT / d y* y- 

i-i Jo 



Nm r 



N 

E 



A 



kp 1 /* T\4 Ap r 



4p r ,N ~ 4p rjA r 



(57) 
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The angular integral An can be solved by applying the following recursion relation: 

A N -i 



An — An -i 



>sm 2N 3 6 cos< 



yielding 



A 



A N -i 



An-2 



N 



2(N- 1) 



An-3 



2(N-l) 2 2 {N -l)(N -2) 2 3 {N -l)(N -2)(N -3) 
Accordingly, by using A 2 — 1/2 one obtains: 

1 



.4 



JV 



2 N - 1 (N-1)\ 



The Eq. ([63]) then reduces to: 



2c 



3/2' 



JV 



1 1 



N\ 2 N - 1 (N~l)l J 



dr r 



2JV-1 



x IV- 



N 



exp ■ 



1 1/3 B 



hp 1 /* T v 4 4 P 



(67) 



(68) 



(69) 



(70) 



(r 2 + Nm c 



Appendix D 

Performing the derivatives, the system in Eq. (|44|) reads 



1 



(71) 



By subtracting the first equation from the second equation, one easily obtains Eq. (I46p . 
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